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Dislocations in soft condensed matter systems such as lamellar systems of poly- 
mers, liquid crystals and ternary mixtures of oil, water and surfactant (amphiphilic 



systems) are described in the framework of continuum elastic theory. These systems 



O ■ 

are the subject of studies of physics, chemistry and biology. They also find applica- 



tions in the industry. Here we will discuss in detail the influence of dislocations on 
the bulk and surface properties of these lamellar phases. Especially the latter prop- 
erties have only been recently studied in detail. We will present the experimental 
evidence of the existence of screw and edge dislocations in the systems and study 
their static properties such as: energy, line tension and core structure. Next we will 
show how does the surface influence the equilibrium position of dislocations in the 
system. We will give the theoretical predictions and present the experimental results 
on thin copolymer films, free standing films of liquid crystals and smectic droplets 
shapes. In semi-infinite lamellar systems characterized by small surface tension the 
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dislocation is stabilized at a finite distance, h eq , from the surface, due to the surface 
bending elastic constant, K s (for zero surface tension h eq K s /2K, where K is the 
bulk bending elastic constant). For large surface tension the edge dislocations are 
strongly repelled by the surface and the equilibrium location for finite symmetric 
systems such as free standing liquid crystal films shifts towards the center of the 
system. The surface is deformed by dislocations. These deformations are known as 
edge profiles. They will be discussed for finite systems with small and large surface 
tension. Surface deformations induce elastic interactions between edge dislocations, 
which decay exponentially with distance with decay length proportional to VD, 
where D is the size of the system normal to lamellas. Two screw dislocations in 
finite system interact with the logarithmic potential, which is proportional to the 
surface tension and inversly proportional to D. The surface induced elastic inter- 
actions will be compared to, well known, bulk deformation induced interactions. 
A new phenonenon discussed in our paper is the fluctuations induced interactions 
between edge dislocations, which follows from the Helfrich mechanism for flexible 
objects. At suitable conditions edge dislocations can undergo an unbinding transi- 
tion. Also a single dislocation loop can undergo an unbinding transition. We will 
calculate the properties of the loop inside finite system and discuss in particular 
the unbinding transition in freely suspended smectic films. We shall also compute 
the equilibrium size of the loop contained between two hard walls. Finally we will 
discuss the dynamical bulk properties of dislocations such as: mobility (climb and 
glide), permeation, and helical instability of screw dislocations. Lubrication theory 
will also be discussed. 

PACS numbers: 61.30.Jf, 61.72.Lk, 61.41. +e, 61.30.Cz 
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I. Introduction 

Before we enter the main subject of this review, which is the edge and screw dis- 
locations in soft condensed matter systems such as smectic liquid crystals, lamellar 
phases of diblock copolymers and lamellar phases in ternary mixtures of oil, water 
and surfactants (amphiphilic system) we would like to put the subject of disloca- 
tions in a broader perspective. The phenomena of upheavals of the mountains, 
slipping of glaciers from the high mountains and easy deformations of metals are 
all controlled by the defects known as dislocations. In the defects free samples of 
solids the shear displacement requires (as theoretically calculated) shears 5 order 
of magnitude larger than is routinely applied in typical deformation processes (for 
an Introduction to dislocations in solids see 1 ' 2 ). This fact aroused the curiosity 
of scientists and in 1936 G.I. Taylor, E.Orowan and M.Polanyi resolved the prob- 
lem of deformations by postulating that deformations of solids are controlled by 
crystal imperfections (defects) called dislocations. It took almost 20 years before 
scientists observed dislocations in solid. One of the first observations were done 
in 1953 by Hedges and Mitchell and later in 1956 by Hirsch, Home and Whelan 
and independently by Bollmann. The two techniques used in the experiments were: 
the decoration technique with optical observation and the transmission electron mi- 
croscopy. In the decoration technique the guest atoms in a solid locate themselves 
along the dislocation line providing contrast for optical observations 3 . In the trans- 
mission electron microscopy (TEM) the diffraction and interference of electrons 4 
allows to see the dislocations. Strong scattering of electrons limits the application 
of TEM to rather thin samples. 

Dislocations also appear in liquid partially ordered lamellar phases(Figs.l-3). 
The molecules in lamellar phases are arranged in liquid like layers parallel to each 
other. Although the symmetry of these phases is the same in amphiphilic systems 5 , 

3 



diblock copolymers 6 and liquid crystals 7 , their structure parameters can be very 
different. In liquid crystals the period of the lamellar structure, d, is set by the 
size of the molecule which is roughly 30A; in amphiphilic system the amount of 
surfactant, water and/or oil allows to vary the lamellar size from ~20A(for large 
concentration of surfactant) to ~1000A(for very small concentration of surfactant) 8 . 
In block copolymers the period is determined by the size of the polymer molecule. 
In the strong segregation limit when the polymer molecules are strongly stretched 
in the lamellas the period is given by M 2 / 3 /, where M is the number of statistical 
segments of a polymer and / is the length of a segment, thus the period can be 
as large as few hundred A. Also the elastic constants characterizing the lamellar 
phases can be very different for these systems and as we will show the behavior of 
dislocations is sensitive to the value of elastic parameters characterizing a system. 
In the following we shall restrict our discussion to uniaxial lamellar phases, i.e. 
smectic A liquid crystals and L a phases of lyotropics. Two types of dislocations, 
namely screw and edge dislocations, are shown schematically in Fig. 4. 

The bulk properties of dislocations have been very thoroughly studied 9-11 since 
1972, when de Gennes presented the equation for elastic deformations induced by a 
single edge dislocation in the infinite lamellar system 7 ' 12 . The first observation of el- 
ementary edge dislocation has been done by Meyer et al in 1978 13 , although in earlier 
experiments the dislocations of large Burgers vector have already been seen 14-18 . 
The influence of surfaces on dislocations has been first studied by Pershan 19 . He 
found that dislocations are repelled from solid surfaces and attracted by a free sur- 
face as in solids. The latter conclusion is erroneous as shown by recent calculations 20 
and experiments 21 ' 22 . Indeed, surface tension effects cannot be neglected at the air 
- smectic interface, contrary to Pershan's assumption. For instance, recent ob- 
servations by TEM (transmission electron microscopy) 22 and AFM (atomic force 
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microscopy) 21 ' 22 of ultra-thin films of A-B diblock copolymers deposited on a solid 
substrate showed that elementary edge dislocations stabilize within the film and 
are repelled from the free surface. Moreover Pershan analysis failed to explain the 
shape of the surface deformations induced by dislocations. The repulsion of the 
dislocation from the surface is also observed in free standing films 23 and in the 
"gouttes a gradins" of Grandjean 24 . The repulsion from the free surface of the 
giant dislocations separating two terraces of a "goutte a gradins" was first pointed 
out by C.Williams 25 in his thesis in 1976. Williams also emphasized that surface 
tension was responsible for this repulsion. 

One of the aims of this review is thus the discussion of the recent progress 
in the dislocation theory for systems bounded by surfaces. The list of problems 
includes: surface deformations (edge profiles), equilibrium location of dislocations 
near surfaces and surface induced interactions between edge and screw dislocations. 
We shall also discuss some of their bulk, static properties such as the core structure 
and the formation of the giant dislocations (which do not exist in solids) and some 
consequences of their flexibility on their interactions. Finally we will analyse the 
dynamic properties of dislocations. These analysis will include the mobility, their 
instability under mechanical stress and their role in microplasticity and rheology. 

The plan of the paper is as follows: In section II we present the experimental 
evidences for the existence of dislocations in lamellar phases. In section III we 
discuss the elastic theory for lamellar phases and in section IV we briefly recall the 
hydrodynamic equations governing their dynamical properties. In section V the 
distortions induced in the bulk by edge and screw dislocations are calculated. Here 
we also calculate the energy and line tension of dislocations. In section VI we find 
the equilibrium location of dislocation in finite system bounded by surfaces. Here 
we discuss both the theory and experiments. The edge profiles are calculated in 
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section VII and elastic interactions are calculated in section VIII. Sections IX, X 
and XI are devoted to the flexibility of dislocation line and unbinding transition in 
the system of dislocations. Finally in sections XII, XIII and XIV the dynamical 
properties of dislocations are discussed. Conclusions are contained in section XV. 
II. Experimental observations of dislocations 

One of the main problem in optical detection of defects is the sufficient optical 
contrast against the lamellar background. In 1978 Meyer et al. 13 observed an array 
of elementary edge dislocations in a wedge shaped sample of smectic liquid crystal, 
contained between two glass plates. In order to increase the optical contrast they 
used the phenomenon of phase transition between smectic-A, where molecules are 
roughly perpendicular to layers, and smectic-C, where liquid crystalline molecules 
are tilted with respect to layers. Tilting is accompanied by the decrease of the layer 
thickness and thus can relieve stress in a wedge shaped sample. Because of the 
wedged shape the thickness of the sample changes and the stress induces formation 
of dislocations in a sample as shown schematically in Fig. 5. For temperatures 
close to the smectic-A smectic-C transition temperature the regions close to the 
dislocations are tilted. Thus under the polarizing microscope one can get dark 
regions separated by bright regions, depending on the tilt and the angle between 
the crossed polarizers. The distance between the dislocations in the array, / is 
given by I = d/a, where a is a small angle between the glass plates. The accurate 
measurements of the angle were done by the Fabry-Perot interferometer. In this 
method one cannot tell wheather dislocations are located at the surface or in the 
bulk, but the fact that only dislocations satisfying the geometrical constraint given 
by / remain after annealing suggests that they are free to move and are not attached 
to the solid surface. This is confirmed by creep experiments discussed in section 
XIII. 
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Another optical method of observation based on fluorescence has been applied 
to phospholipid lamellar phase (example of self assembling amphiphilic system) 
with alternate layers of water and fluid-like bilayers of phospholipid 26 . The system 
is placed in the same wedge shaped container and the phospholipid is doped with 
fluorescent, lipid analogue molecule. The sample is annealed for some time (several 
weeks) and next excited with a laser. The fluorescent intensity is monitored as a 
function of position as shown in Fig. 6. The visible steps are clear signs of parallel, 
elementary edge dislocations. The long annealing time that is necessary to reach 
equilibrium in this system is related to permeation (flow of molecules across lamellas, 
Sections XII, XIII) . Permeation is much more difficult in lyotropic systems like this 
one than in thermotropic systems studied by Meyer et al 13 , where the annealing 
times were much shorter. 

The exact location of dislocation in a sample can be determined by TEM 22 . 
Here we briefly describe this experiment 21 ' 22 as shown schematically in Fig. 7. The 
ultra-thin films (thickness D ~ 1000A) of lamellar phase of AB diblock copolymers 
(period d ps 315A) are deposited on solid substrate. The films are first prepared 
in the disordered solid state and heated above their glass transition but below 
the order-disorder transition. The lamellar (smectic) ordering appears with layers 
parallel to the substrate. The free surface of the film is not flat since circular steps 
of micrometer size are formed, separated by equally large holes. The process of the 
formation of the holes and steps is associated with the nucleation of dislocations 
inside the film and is due to the incompatibility of the initial thickness of the 
disordered film, D, with the smectic period of the final lamellar phase i.e. Nd < 
D < (N + l)d. Because N + 1 is the number of layers below the step and N is 
the number of layers below the hole the height of the step is equal to the lamellar 
period, d. Indeed applying TEM one can observe the dislocation itself, while by 
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applying AFM one can measure the shape (size, height) of the step and the width 
of the profile at the edge of steps (edge profiles). The slices (cross section shown in 
Fig. 7) of thickness of 500 A are prepared with layers parallel to the electron beam 22 . 
To have a contrast between the PS and PBMA domains PS are stained by treating 
it with RuC>4 vapor. The electrons are diffracted on the layers and the micrographs 
shows the positions of layers and in particular of defects. This is a direct way 
of seeing dislocations. The slices have to be thin, because of strong interactions 
between electrons and matter. The originally liquid samples are of course vitrified 
before observation. Since the polymers are very viscous fluids, the structure remains 
intact in the process of vitrification. 

Edge dislocations are also visible in free standing films 23 . They are usually 
obtained by smearing a smectic phase with a wiper across a fixed aperture made 
in a solid substrate. This technique used by most researchers and invented by 
Friedel 70 years ago 27 has been considerably improved by P.Pierahski 23 . He used a 
rectangular frame of variable area, allowing precise control (up to one layer) of the 
film thickness in the process of preparation. Fig. 8a shows a schematic structure of 
the film near the edge of the aperture. As we see the dislocations are located in the 
middle of the smectic meniscus. Fig. 8b shows a typical arch - texture obtained when 
the film is quickly stretched. Here the film is divided into areas of different thickness 
which are separated by lines of circular shape. In general the thickness variation 
between two adjacent regions is one layer. Some lines are linked by knots too. 
Pierahski 23 showed that these lines are bulk dislocations located in the middle of 
the film and suggested that there exists some dust particles (not always visible under 
the microscope) at each junction (the knots) between the dislocations (Fig. 8c). 

Edge dislocations also play a crucial role in the problem of the equilibrium shape 
of a smectic A droplet deposited on a glass plate with strong homeotropic anchoring 
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(molecules arranged perpendicularly to the substrate) 28 ' 29 . Droplets shapes can be 
measured very precisely (±50 A) by Michelson interferometry. This experiment was 
performed with materials that have either a nematic-smectic A phase transition 
(80CB, 40.8) or an isotropic-smectic A phase transition. It revealed that small 
droplets (usually less than 200 micrometers in diameter) have a single facet parallel 
to the substrate, while larger droplets are shaped like a spherical cap. Fig. 9 shows 
typical profiles of a small 80CB droplet, when the temperature is decreased below 
Tna (nematic-smectic phase transtion temperature). Because this transition is sec- 
ond order, the size of the facet vanishes at Tna and increases like (Tna — T) 45 
when the temperature is decreased below T N a- In 10OCB where the first order 
smectic-A isotropic phase transition occurs, the facet has a finite size at the tran- 
sition. If the radius of the droplet is larger than 200 micrometers, the top facet 
disappears and the droplet becomes spherical as in the previous case. It means that 
steps on the free surface associated with facets (Fig. 10a) are replaced by the disloca- 
tions distributed throughout the interior of the droplet (Fig. 10b). This is possible 
since dislocations are repelled by both, the solid substrate and the free surface. 
In very large droplets, dislocations can group together and form giant dislocations 
separated by terraces parallel to the glass plate ("gouttes a gradins" Fig. 11) 

The shapes of the droplets has been explained by the theory, where interactions 
between the steps and between the dislocations, together with the surface and finite 
size effects have been properly taken into account 30 . These finite size efects are also 
important to explain the appearances in some materials (40.8) of a secondary facet 
around the primary one (Fig. 12). This facet occurs at low temperatures and in 
large enough droplets (diameter more than 100 micrometers). This means that 
very special conditions must be fulfilled simultaneously in order to observe it. This 
unexpected facet (there is no terrace periodicity along the new facet direction) 
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results from the elastic attraction between like steps, balanced by entropic repulsion. 
We shall see in section VIII that there can be an attraction between like dislocations 
(of the same sign) in films, similar to the attraction between steps as observed in 
this experiment. 

The screw dislocations have been also observed in self assembled amphiphilic 
systems in TEM 16 ' 31 ' 32 experiments (Fig. 13). The density of screw dislocations 
can be as large as 10 8 per cm 2 in many lyotropic systems. High density suggests 
that their energy is rather small, a fact also established in the theory (Section V). 
III. Elastic theory of lamellar phases 

The lamellar phase consists of parallel, equidistant two dimensional liquid lay- 
ers. The average distance between the layers is denoted d. For convenience we 
assume that unperturbed layers are perpendicular to the z-axis. The layers de- 
form easily and their deformations are conveniently described by two unit vectors, 
one normal to the layer, m(r) and one describing the average orientation of the 
molecules n(r) (called in liquid crystals the director), and a scalar quantity, which 
measures the distance between the layers along the normal. We further assume 
that other quantities like the nematic order parameter and the density (for liquid 
crystals) or surfactant and water concentrations (for amphiphilic system) adjust to 
the layer deformations. In the first approximation, density changes do not gener- 
ate new terms in the distortion energy, but only renormalize the lamellar elastic 
constants 7 ' 33 . Using the two aforementioned quantities we can write the distortion 
energy density in the simplest form invariant with respect to global rotations: 

' (3-1) 
+ -K 3 \n(r) x (V x n(r))| 2 + ifl^n - m) 2 , 

where do is the unperturbed layer spacing, B is the Young modulus of the layers 
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associated with the variations of the thickness and Bi ~ B is the elastic modulus 
associated with the deviations of the average orientation of molecules from the 
normal to the layers 34 . The divergence and the curl operators in Eq.(3.1) are taken 
with respect to variables in the distorted state and the point r is uniquely related 
to some other point ro on the unperturbed layer. We neglect the higher order 
terms in d — do and in the derivatives of n. The last three terms in this distortion 
energy have their analogs in nematic liquid crystals. K 1 , K 2 and K 3 are the elastic 
constants for splay, twist and bend deformations of n, respectively. As one can see 
the splay of n corresponds to the bend of the layers while the bend of n corresponds 
to the splay of the layers. From now on we shall use the notion bend and splay in 
connection to layers. For the last term in Eq(3.1) one finds that (n — m) 2 ~ X/R 
where A = ^jKi/B is roughly proportional for most of the lamellar systems to do 
and R is a typical radius of curvature for the director field. Consequently for R 3> A 
we can set n = rh. We shall keep in mind however that this assumption is not valid 
close to the dislocation core. Nonetheless from now on we neglect the last term in 



Since at this point we are interested in the local distortion energy density in 
the lamellar phase we can neglect dislocations, which means that the total number 
of layers crossed along any path going from some point A to another point B is 
constant i.e. for a closed path 



where <i, in general, depends on the position r. Condition(3.3) eliminates the twist 
term from the distortion energy, but not the splay term (K3). Here we only note 



Eq(3.1). 




(3.2) 



This is equivalent to the condition 7 : 



vxH-o. 



(3.3) 



11 



that splay of layers does not affect the long wavelength properties of the system. 
Combining Eq(3.1) and Eq(3.3) we get 



+ Ki|Vn(r)| 2 + K 3 



d 



Vd 



n(r) 



d 



Vd 



(3.4) 



It is important to notice that a dislocation should not affect the form of the distor- 
tion energy density, which is a local quantity, whereas dislocation is described by 
the global condition i.e. nonvanishing of (3.2). One can also verify this by direct 
calculations 7 . For the edge dislocation along y direction n has only x, z components 
and does not depend on y and V x n/d is along the y direction. Thus the twist 
term drops out in this case. 

In lamellar phases, n and d are not independent and can be expressed in terms 
of the vertical displacement field u(tq). Although u is a multivalued function in the 
presence of dislocation, we can still use it away from the dislocation. The vector 
normal to the layer at point r = r(ro) is simply given by: 



n = 



du _ du y 

dx ' dy ' 



(3.5) 



y^l + ivf u\ 2 

where r = (^cyo^o) * s the coordination point in the unperturbed system. The 
distance between layers measured along n is 



d = d z (hz), 



(3.6) 



where d z is the distance between layers measured along the z-axis. We find, to the 
lowest order in the derivatives of u, that 

du 



d z = d 1 + 



dz 



(3.7) 



which together with Eqs(3.5-6) gives d in terms of tt, 



d = do 



1 + |vfw| 2 



(3.8) 
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In order to represent the distortion energy given by Eq(3.4) in terms of u, we have 
to transform the point r on a perturbed layer to the point r on the unperturbed 
layer and the corresponding nabla operator as follows: 



X = Xq, 



y = 2/0, 



z = z + u(x ,yo,zo), 



and 



d_ 

dx 

d 



ddu 
O 

dx Q + du_ j dz ' 



d 



du 
dyo 







dy dy ( 1+ du_\ dz ' 

V dzo J 

1 d 



d_ 

Yz 



(3.9) 



(3.10) 



(3.11) 



Combining these equations we recover the condition given by Eq(3.3) for the con- 
stant number of layers. Finally we obtain the following expression for the distortion 
energy density 35 : 



fb 



3" 



d — do 
d 



+ K, 



7(0), 



+ K 3 



d 



__ (3-12) 

where nj_ = — Vj°^w/ \J 1 + | Vj°^w| 2 and Vj°^ = (d/dx , d/dy ) is the two dimen- 
sional nabla operator. As noted before, although the splay term associated with K3 
gives a coupling between compression and undulation of layers, it does not affect the 
long-wavelength properties of the system. Also it can be neglected in comparison 
to the compression term. Ignoring the splay term and anharmonic contribution to 
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the bend term, we find an approximate form of ft,: 
/ / 



fb 



B 



1 + 



du(r ) 

dzn 



II +K 1 (Af«(r )) 



V 



l + |vf u(r )\ 2 



\ 



I 



(3.13) 



Here is the two-dimensional Laplacian with respect to x ,y . The distortion 
energy density given by Eq(3.13) is invariant with respect to rotations. For example, 
by rotating the system around y-axis by an angle 6>, we find: 



u = z 



1 



cos 9 



— 1 ) — xq tan 6. 



(3.14) 



For such spurious deformation we find from Eq(3.13) that fb = as it should. Ex- 
panding Eq(3.13) in iV^wl 2 we recover the de Gennes distortion energy density 7 ' 36 



-i|vf«(ro)| 2 ) 2 + ^(Ai°)«(ro)) 2 j, 



(3.15) 



in the lowest order of the expansion, losing however the rotational invariance of fb- 
Alternatively, the distortion energy density can be expressed in terms of Vw. We 
find 



f = - f B I 1 

' 2 1 I v /l-2(du/dz-\X7u(r)\ 2 /2) 



1 +Ki|Vn(r) 



(3.16) 



where 



n(r) = 



{—du/dx, —du/dy, 1 — du/dz) 



(3.17) 



'I- 2 (du/dz - |Vu(r)| 2 /2) 

and we have neglected the splay term. One easily checks that = when we rotate 
the system, i.e. 

u = z(l - cos 6) - xsin6. (3.18) 

The transformation properties of u result from the fact that it is defined as the 
vertical displacement of the layer from its rest position. In the lowest order expansion 
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in du/dz — \Vu 



/2 we recover the Grinstein and Pelcovits expression 



'•-if 8 



( 



du(r) 1 
~dz 2 



V«(r)| 2 +K 1 (A ±W (r)) 2 



(3.19) 



If we neglect the anharmonic terms we find the harmonic approximation of the 
distortion energy: 



Here we define K=K 1 . Once again we note that Eq(3.20) is valid far away from 
the dislocation core. This form of the bulk distortion energy will be used in our 
description of deformations induced by dislocations in a lamellar (smectic) system. 
At the harmonic level of approximation the form of the distortion energy is the same 
in coordinates of both, the deformed and undeformed state (Eq(3.15), Eq.(3.19)). 
The role of the anharmonic terms in lamellar phases for long wavelength distortions 
has been discussed in detail in Refs(7,37). In general its role for distortions induced 
by dislocations is not known. 

IV. Hydrodynamics of smectics: elastic and viscous stress tensor 
The motion of dislocations can be described in the framework of the hydrody- 
namical equations for the lamellar phases. Here we recall the basic ingredients of 
the theory 7 ' 33 . Let v be the average velocity of the molecules. If the velocity of the 
dislocations is much smaller than the velocity of the first sound, then the medium 
can be regarded as incompressible i.e. 




(3.20) 



Vv = 



(4.1) 



The momentum conservation equation for the system is as follows. 



Dw 



= Va + F, 



(4.2) 
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where p is the density, a is the stress tensor and F the bulk forces (gravitation for 
example). The stress tensor can be decomposed into three terms: 

a = -PI + a E + (7 V (4.3) 

Here P is the hydrostatic pressure (given by condition (4.1)), cr E is the elastic 
stress tensor and a v is the viscous stress tensor. The former is related to the layer 
displacement u and can be obtained from the following equation relating the change 
of the free energy fb to the displacement: 



du 



Since u is the z-component of the displacement vector, the only nonvanishing compo- 
nent of this tensor have (z,j) indexes. Expanding (3.20) and equating the expansion 
to the right hand side of (4.4) we find 

a E = _^J± _ _d__dfb_ , 45 ^ 
ZJ duj dn du,ij ' 

Here uj and the first and the second derivatives of u with respect to the 

components of r. More explicitely we find: 

= ~K§- x (Aiu) , (4.7) 

< = ~ K §y ( A ^) ■ ( 4 - 8 ) 

All other terms are zero. Note that elastic stress tensor is not symmetrical due to 
the curvature elasticity and the associated surface torque Cn . The nonvanishing 
components of the tensor C are 

C xy = (4.9) 
ua ,yy 
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and 



C„ = n (4.10) 



The elastic stress tensor satisfies the torque balance equation 

Cij,j - e ljk af k = (4.11) 

automatically (we assume as stated in section III that n is normal to layers). 

In order to calculate the expression for the viscous stress tensor we have to 
find the thermodynamic fluxes and forces associated with the entropy production. 
One finds for the isothermal processes the following expression for the change of the 
entropy s in time: 



where 



G = <r? jd (4.13) 



is the elastic force normal to layers and 

*-H§ + £)- (4 - 14) 

Now from the Onsager relations between thermodynamic forces (G, ajj) and fluxes 
(Aij, (Du/Dt — v z )) we find 

ajj =ixi5ijA kk + fi28iz8j z A zz + ^Aij 

(4.15) 

+^4 (SizAzj + S jz Azi) + H^iz^jzAkk + SijAzz) 



and 



^-v z = X P G. (4.16) 



Here jii (i = 1 • • • 5) are independent viscosities. In the following we shall take 
ajj = [iAij for simplicity. Eq.(4.16) is formally equivalent to the Darcy law in 
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porous medium. In our case the lamellar phase plays both the role of the fluid and 
of the porous medium. It can be shown that the permeation coefficient X p , 

Xp ~ ~i^r (4 - 17) 

where D\\ is the diffusion coefficient normal to the layers and v m oi is a molecular 
volume. The permeation coefficient is much smaller in lyotropic (e.g. amphiphilic 
systems) than in thermotropic smectic liquid crystals. This coefficient can be esti- 
mated from the measurements of the edge dislocation mobility as will be shown in 
Section XIII. 

V. Static bulk properties of dislocations 

Let us consider an elementary edge dislocation (of unit Burgers vector) located 
along the y-axis at x = I and z = h. The displacement u is a multivalued function 
in the z = h plane i.e. 

/ \ _ / 0, if x < I; , , 

u( - X ' z) = \sgn(z-h)d/2, ifx>l. {bA) 

This is equivalent to the condition 

^1=1, (5.2) 

where the contour of integration is around the dislocation line. The distortions 
induced by a dislocation are described by the following equation 7 : 

- X 2 A 2 L u = (5.3) 

obtained from the minimization of the total bulk distortion energy i^, 

F h = J drf b . (5.4) 
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Here ft, is the energy density given by Eq.(3.20) and A = ^JK/B. The equilibrium 
solution of Eq.(5.3) satisfying condition (5.1) is as follows 7 : 



Ub(x,z) =sgn(z — h) ^— J dqexp (—Xq 



2i z _ h \\ ex P( i( l( x 



i(q-iO+) 

= -sgn(z — h) 1 + erf 



4 \ \2y/X\z-h\ 

where erf(rj) = J * exp(—u 2 )du is the error function. The solution is valid 

outside the dislocation core of size r c . For dislocations of Burgers vector, b, of 
length nd, for which 

^dl = in, (5.6) 

the solution is simply given by ±nub (here n is a positive integer). Inserting the 
solution (5.5) into Ft, (Eq.(5.4)) and introducing the cutoff r c in (5.4) we find the 
energy of the edge dislocation line per unit length F b /L y : 



Eo = Ec + (5 . 7) 

2r c 

where E c is the core energy. As expected this energy does not depend on h and / 
in the infinite systems, but as we shall see in the next section it does depend on the 
location of dislocation (/, h) in the finite system bounded by surfaces. The energy of 
dislocation per unit length is finite in lamellar phases whereas in solids it diverges 
logarithmically with the size of the system. Also stresses do not vary like 1/r as in 
solids. We find from section IV (Eq(4.6)) and Eq(5.5) that in lamellar phases 

= { *~ l \ d . u exp (-(* - /)V(4Ak -h\)). (5.8) 

o\z — h\y'nX\z — h\ 

Thus af z is large only for 4A|^ — h\ > x 2 . Outside this region the stress vanishes 
very quickly because of the fluidity of layers. Let us now discuss the core structure 
of an edge dislocation. Let us assume that the core is nematic. We also assume 

19 



that the core is anisotropic i.e. along the z direction it has the size 2r c while along 
the x direction it has negligible size. Then 

E c = 2 7 r c (5.9) 

where 7 is roughly proportional to the nematic-smectic surface tension. If we now 
minimize Eq (Eq(5.7) with respect to r c we find 

E = 2^B\^nd (5.10) 

and 

nd BX 

r « = TVT- (6 ' n) 

When 7 is small the core size can be very large as shown in (Fig. 4b). The total 
energy is now proportional to the length of the Burgers vector, nd, of an edge 
dislocation and this explains why dislocations may group together. When their 
Burgers vector is large (giant dislocations with n > 10), the core of dislocation 
can split 17 ' 11 ' 7 into two ±1/2 disclinations a distance nd/2 apart (Fig. 4c) The core 
energy is in this case, 

E c = — In — )+w c , (5.12) 



2 V> , 

where w c is the disclination core energy and tq is the disclination core size. Now 
the core energy is lowered by the gathering of elementary edge dislocations into 
dislocations of large n. The core structure of disclinations in smectics is not known, 
although recently the disclination core structure in nematic liquid crystals has been 
thoroughly studied in the Landau-de Gennes model 38 and in computer simulations 
for hard rod systems 39 . The core structure of this defect depends on the length of 
molecules. For short molecules, the core is biaxial with molecules perpendicular to 
the disclination line, whereas for long molecules they are parallel to this line. The 
core size is proportional to molecular thickness, rather than molecular length. 
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Let us now discuss the screw dislocation properties. In this case (Fig.4d) the 
Burgers vector is parallel to the dislocation line and perpendicular to layers. Recent 
discovery of the TGB (Twist Grain Boundary) phase consisting of lattice of screw 
dislocations 40 ' 41 (analogous to the type II superconductors) brought the renewal 
interests in these defects 42 . The static properties of screw dislocations have been 
intensively studied 42-47 and here we shall briefly present the major results of this 
properties. The distortions induced by screw dislocation of strength n located along 
the z axis (perpendicular to smectic layers) are also described by Eq.(5.3) and the 
condition(5.6). For small distortions we have: 

V±ud\ = ±nd, (5.13) 

Due to the symmetry u does not depend on z and the layer thickness d is equal 
to the undistorted value d . We find in this case the equilibrium solution of the 
form 43 : 

u b (x, y) = ±^ = (arctan (|) - tt + (7r/2)sgn(y)) . (5.14) 

where <p is the polar angle in cylindrical coordinates. The solution is only valid far 
away from the core, i.e. for p ^> r c , where r c is the core size and p is the distance, in 
polar coordinates, form the screw axis (given by z = 0). The total distortion energy 
is the core energy only, since the strain energy is zero. In solids the strain energy 
diverges logarithmically with the size of the system. If we include anharmonic 
terms in the distortion energy (section III, Eq(3.15)) and the splay term, which for 
distances not too far from the core can compete with the compressional term, we 
find the energy per unit length of the screw dislocation 43 ' 46 in the following form: 

E - E I {nd) ' B I {nd) ' K * (5 151 

It should be noted that we have used the solution (Eq.(5.14)) which followed from 
the harmonic approximation (Eq.(5.3)). Since K 3 and Bd 2 are of the same order of 
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magnitude, so both contribute significantly to the distortion energy. Note however 
that this elastic energy is three orders of magnitude smaller then the energy of the 
edge dislocation. The core energy (core size r c ) can be roughly estimated 44 as: 

E c = k B (T NA - T)nr 2 c /v moh (5.16) 

where v mo \ is the molar volume and T^a is the transition temperature from the ne- 
matic liquid to the lamellar phase. The core energy vanishes at the transition which 
means that the formula (5.16) is valid only for the second order phase transition. 
In the case of the first order phase transition we find: 

2 AH(T N a -T) , m m 
E c = nr c + 2r c7 , (5.17) 

J-NA 

where AH is the latent heat (per unit volume) of the transition and 7 is the nematic- 
smectic surface tension. Here it has been assumed that the core is filled with 
undercooled nematic phase. Similarly as we did in the case of edge dislocations 
(Eqs(5. 10-11)) we could minimize Eq with respect to r c . This calculation shows 
that the energy of the screw dislocation increases faster than the length of the 
Burgers vector nd so that elementary screw dislocations are favored. In fact the 
exact core structure of screw dislocation is not known in general, although some 
preliminary studies have been done for smectic liquid crystals 42 ' 46 . These studies 
indicate that the core is nematic (although there is a possibility for the completely 
isotropic core) . Since layers are not well defined inside the core one has to use there 
the Landau-de Gennes description in terms of the order parameters 7 . The smectic 
order parameter e, the nematic order parameter s and the local angle between the 
director field and the z axis, are shown in Fig. 14. The nematic order parameter 
does not change in comparison to its bulk value S}>. The smectic order parameter 
decreases to zero in the core as p n , and approaches the bulk value e& as — 1/p 4 . 
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Within the core the director n tends to lie along the screw axis and $ ~ p close to 
the center of the core. The director n = e p cos sin 1} + sin a b sin 1} + e z cos 1} 
should be identified far from the core with the layer normal. Far from the core one 
finds ctb = 7r/2 and 1} = arctan(n<i/27rp). The size of the core increases with the 
strength of the screw dislocation n. Finally we note that the internal stress inside 
the core (giving rise to radial force) is balanced in the smectic by the non-uniform 
pressure p = p + (nd) 4 B / '128n 4 p 4 , which can be as large as twice the athmospheric 
pressure, po, of a defect free sample. The case of the isotropic screw core structure 
has also been studied 42 ' 46 . These cores are energetically unfavorable, although they 
may be stabilized close to the smectic liquid crystal - isotropic phase transition. 
The core structure for the polymers or amphiphilic systems certainly cannot be 
nematic since the nematic phase is not present in these systems. To our knowledge 
these cores have not been studied in detail (see however Refs.11,44). 

Finally let us mention an important property of screw dislocations, namely the 
fact that although their energy per unit length is small their line tension is large 44 . 
In order to calculate the line tension of the screw dislocations let us assume that 
the dislocation is tilted with respect to the layers, making an angle 9 with the z 
axis Fig. 15. For very small 9 the displacement is given by: 



Ub(x, y, z) = — arctan (5.18) 

2% x 



and the total energy is 



E (9) =E (9 = 0) + ln (5 JQ) 

where R is the size of the system or the distance between the neighbouring screw 
dislocations of opposite signs. The line tension is defined as follows: 

T screw = £ (# = 0) + ^|^ (5.20) 

d9 2 9=0 
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and it follows from (5.19) that T screw 3> E contrary to the case of edge dislocations 
for which T e d ge ~ E . The main difference between the two types of dislocations is 
that edge dislocations located along the y axis do not change their character when 
tilted in the x-y plane with respect to the y-axis, while the screw dislocations along 
the z axis acquire the edge character when tilted with respect to the z axis. This 
result explains why screw dislocations prefer to be perpendicular to smectic layers 
and do not tilt easily (Fig. 16). We shall also see that this property is crucial for 
understanding the helical instability under compression or dilation normal to layers. 
VI. Equilibrium location of dislocations in finite systems 
Experimental systems are always bounded by surfaces. In Meyer et al 13 ex- 
periment two glass plates bound the lamellar phase (Fig. 5). Thin freely suspended 
smectic liquid crystal films 23 (Fig. 8) are attached to the aperture and bounded 
from above and below by a smectic free surface (smectic-air interface). Smectic 
droplets 29 ' 30 deposited on a solid substrate (Fig. 10) are bounded by the solid sur- 
face and the free surface. The diblock copolymer lamellar film deposited on a solid 
substrate shown in Fig. 7 is also bounded by a free surface. In the finite system 
dislocations are influenced by the bounding surfaces. The boundary effects can be 
incorporated in the form of the surface distortion energy. It should be noted that 
contrary to the solids, the surface energy of distortion in lamellar liquid phases 
is comparable to the bulk energy. In solids the energy scales for the surface and 
the bulk are well separated with surface energy being much smaller than the bulk 
energy. 

The influence of surfaces on dislocations has been studied by Pershan 19 . How- 
ever his analysis was not sufficient to explain interactions with a free surface 25 . 
The problem has been first analyzed theoretically by Lejcek and Oswald 20 in the 
framework of the image dislocation approach. In the semi-infinite system bounded 
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by the free surface at z = the total energy contains the bulk part (5.4) and the 
surface part 20 : 

F s = lJ dr ± <y\Vu s \ 2 , (6.1) 

where 7 is the smectic-air surface tension and u s = u(rj_,z = 0) is the surface 
displacement. Minimizing Fb + F s gives Eq.(5.3) and a boundary condition at the 
free surface: 

B—= 7 A ± u s . (6.2) 

Let us now consider a simple case of a dislocation of Burgers vector b = nd located 
at z = h. It can be shown that both equation (5.3) and (6.2) are satisfied by the 
solution (5.5) if we introduce a symmetric image dislocation (located at z = —h) of 
Burgers vector: 

bi = b (6.3) 

and take the linear combination of the bulk displacement field (ub Eq(5.5)) for both 
dislocations u = nubih) + niUb(-h). Note that bi = nid does not have to be a 
multiple of d since it is a virtual dislocation and consequently Hi in general is not 
an integer. We shall see in section VIII that two dislocations repell each other when 
have the same sign and attract each other for opposite signs. Thus it follows that 
the surface attracts dislocation for 7 < v 'KB and repels it in the opposite case. 

These calculations can be generalized to more complicated situations, when for 
example we have two surfaces 20 ' 30 , but in this case we have to introduce an infinite 
number of image dislocations, which is not very convenient. For this reason we shall 
adopt another approach to the problem based on de Gennes method of calculation 7 . 
We also note that the surface energy considered by Lejcek and Oswald is incomplete 
since it does not contain the surface curvature term which influences very strongly 
the equilibrium location of dislocations in the case of small surface tension 48 ' 49 . The 
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surface distortion energy with the curvature term has the following form 50 

F s = lJ dr± {l\Vu s \ 2 + K s (A ± u s ) 2 ) . (6.4) 

Here the new quantity is the surface elastic constant, K s . It should be identified 
with the bending elastic constant of the last lamellar layer (at the surface). In 
general the properties of the surface layer are different from the bulk layers and the 
value of the surface constant can be different from the bulk value. In particular, 
if the surface layer is the same as the bulk one, they are related by the formula: 
K s = Kd, where d is the layer spacing. The fact that the surface bending elastic 
constant must be included follows from the comparison of the distortion energy in 
the discrete 51 and continuous 50 ' 52 representation. In the discrete representation the 
layers have indexes i = 0, 1, 2... and u(x, z) = Ui{x). Then the last term in Eq(6.4) 
appears naturally with K s = Kd. In the discrete representations we sum over i 
instead of integrating over z and consequently each layer, including the surface one, 
has the bending elastic constant equal Kd. 

Let us consider the experimental system shown in Fig. 7 with the substrate 
located at z = D and an elementary edge dislocation (n = 1) located at z = h and 
x = 1. The condition for the edge dislocation is given by Eq.(5.1). Now minimizing 
the total energy F = Fb + F s (Eq(5.4) and (6.4)) we find Eq.(5.3) plus the boundary 
condition at the free surface (z = 0): 

-B—- 1 A ± u s + K s Alu s = (6.5) 

and the boundary condition at a solid substrate, 

u(x,z = D) = 0. (6.6) 

The equilibrium solution satisfying Eqs(5.1,5.3,6.5,6.6) is given by u eq {x, z), 

u eq (x, z) = Ub(x, z) + u p (x, z) (6.7) 
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where the singular bulk part of the solution, it&, is given by Eq.(5.5) and the partic- 
ular non-singular part of the solution, necessary to satisfy the boundary conditions 
(6.5-6), is 48 

«p(z, *) = ^ / dq eW {lQ ^ ~ l)) (exp (- + exp (Xq 2 z)f 2 (q)) . (6.8) 



Here 



. , v ai(g)e-^ ah + ai(g)e^ a ( fe - 21? ) _% 
a 2 (q) + ai(q)e 2A i u 

a^e-^^D) _ aMe x q *( h -2D) 

hW = t-—. / \ _2A„2 D (b.iu; 

a 2 (q) + ai{q)e 2A <? u 



and 



a 1 (q) = l-a s -X 2 s q 2 , (6.11) 

a 2 {q) = l + a s + X 2 s q 2 . (6.12) 

Here a s = ^/y/KB is the dimensionless constant, whereas A = yjK/B and A s = 
\J K s /\/kTB are two microscopic lengths. The total distortion energy F = F b + F s 
can be written in the following convenient form: 

F(h) =E + ^B J dx dUp ^ Z ^ (u b (x, z = h~)- u b (x, z = h+)) . (6.13) 

Inserting u b and u p into Eq(6.13) gives F(h) as a function of h in the following 
form: 

r , VKBd 2 f -a 1 (q)e- 2 ^ 2h + a 2 (q)e 2X i 2 ( h - D ^ -2 ai (q)e- 2 ^ 2D 

F(h) = E H / dq — — - 2Xa 2 D • 

8tt J a 2 (q) + a^qje 2Xq D 

(6.14) 

As one notices the equilibrium location of the dislocation is still to be determined 
The equilibrium location of the dislocation, h eq , corresponds to the minimum of 
F(h). 
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First let us consider the case of the infinite thickness of the film D = oo. The energy 
as a function of h is shown in Fig. 17. For a s = we find h eq = 1.08K S /(2K) and 
for a s < 0.2 we get roughly h eq > (K S /(2K))((1 + ct 8 )/(l — a s )). For larger values 
of a s , we find that the following inequality holds: 

1 + a s 2h eq K 1 

°- > > . 6.16 

1 — a s K s 1 — a s 

For small surface tension the dislocation is stabilized at a finite distance due to 
the surface bending elastic constant K s . Finally for a s — > 1 we have h eq — > oo, 
which means that the dislocation is no longer stabilized at a finite distance from the 
interface. This result is the same as obtained within the image dislocation approach 
of Lejcek and Oswald 20 . For large surface tension a s ^> 1 and finite but large D the 
dislocation is stabilized roughly in the middle of the system at h eq ks _D/2 19 > 20 > 22 > 48 . 
For small surface tension and finite D the location of dislocation in a system depends 
strongly on D. (Fig. 18) For small K s the dislocation is close to the interface, while 
for larger K s it moves towards the center of the system. For thin films there is a 
strong deviations from the linear dependence of h eq on K s . This is of course due 
to the strong repulsion of the dislocation from the solid substrate, which cannot be 
compensated by any finite value of K s . Finally dislocations in free standing film 
shown in Fig. 8a, are located exactly in the middle of the film, since 7 3> \JKB in 
this case 53 and the film is symmetric. 

Usually the surface tension for the air-lamellar interface is about 7 = 30 dyn/cm 
(this value characteristic for air-organic liquid interface (free surface)) and in gen- 
eral larger by an order of magnitude than \/KB. However, there are also many 



systems characterized by a small value of the surface tension (smaller than yKB). 
They include nematic-smectic interface (for liquid crystals), microemulsion- lamel- 
lar interface in amphiphilic systems and AB diblock copolymer lamellar- A (or B) 
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homopolymer melt interface. Our results can be used to describe dislocations in all 
these systems. 

Here we may also ask the question concerning the average distance, /, between 
dislocations in a smectic meniscus (Fig. 8a). Neglecting the elastic interactions we 
find 

I = Tggg(gZ|) , (6.17) 

Iaw - Ilw - F{h)/d 

for large film thickness and far from the aperture (x > £o)- Here "fAW is the air- 
aperture surface energy(per unit area), ^lw is the lamellar-aperture surface energy 
and F(h) is the energy of dislocation (per unit length) and h = D /2 is half of the 
thickness of the meniscus close to the aperture. £o = \flTpg is the capillary length, 
p is the density and g is the gravitational accelaration. This length can be very 
large (of the order of millimeters). Close to the aperture the dislocations are very 
close, interact strongly and the formula (6.17) may not be valid. That is why the 
smectic meniscus offers a very good opportunity to study the competition between 
gravitational forces and elastic interactions between dislocations. The system can 
also develop sharp step defects at the surface close to the aperture. These steps do 
not have to be related to dislocations inside the meniscus (see discussion on steps 
in Sec. II and droplets in Fig.lOab). The meniscus also plays an important role for 
explaining the tendency of the film to get thinner. Indeed, the pressure inside the 
film P is less than the athmospheric pressure, Pq, because the meniscus is curved. 
The equilibrium condition at the surface gives P — af z = P and thus the layers 
are compressed. Consequently the film tends to get thinner. Let us assume that we 
have nucleated a dislocation loop, say by suddenly pulling a film. Then we find that 
the dislocation loop of size R and Burgers vector of length nd will spontaneously 
grow for R> R c = F(D/2)/((Pq — P)nd). Here D is the thickness of the film. 
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Finally let us discuss the clustering of like dislocations (Fig.l9ab) in a wedge 
shaped samples 13 ' 26 ' 54 ' 55 . Calculation of the energy in the case shown in (hatched 
region in Fig. 19a) gives, for dislocations of Burgers vector nd: 

£ = 2£o(n) + ^! (6.18) 

where a is the wedge angle, and E (n) is given by Eq(5.7) and D is the average 
thickness of the sample over a distance between the neighbouring dislocations. For 
the case shown in Fig. (19b) when dislocations are grouped in pairs of the doubled 
Burgers vector 2nd we find: 

E = E (2n) + (6.19) 

where E (2n) is given by (5.7). If 2E Q (n) > E{2n) then we find that dislocations 
form pairs for the angle, a, larger then the critical angle ct c , 

B{ndf 

ac ~ 4D(2E (n) - E (2n)). ^ 

Experimentally 54 one finds a c ~ 10 -3 when D = 100/Um and n = 1. Then from 
(6.20) we find 2E (n) - E (2n) ~ 10~ S K, which means that 2Eo(n) ~ Eo(2n). Al- 
though they are not equal exactly as predicted by (5.10) they are very close, indeed. 
Certainly the small differences can come from the subtle changes in the core struc- 
ture and/or anharmonic terms. Similar studies have been done for ferrosmectics 56 . 

VII. Edge profiles 
Inserting h eq into u eq we get equilibrium distortions in the film; the edge pro- 
file at the interface is given by u eq {x,z = 0) (see Fig. 7). The influence of the 
surface tension on the width of the profile is shown in Fig. 20. Large surface ten- 
sion broadens the edge profile 22, and for a symmetric film of size D the width of 
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the profile, Ax scales as 22 y^yDjB. Theoretical analysis 22 ' 48 explained qualita- 
tively and quantitatively the results of experiments 21 ' 22 performed on thin diblock 
copolymer films deposited on solid substrate. The comparison between the theory 
and experiment 23 is shown in Fig. 21. The values of parameters obtained were as 
follows: 7 = 23.5 dyn/cm 2 (see also Ref.57), A = y/K/B « d = 315A, yfKB = 0.9 
dyn/cm, 'y/y/KB 26. Also the scaling law for the width i.e. Ax ~ \[D has 
been verified experimentally. The same analysis has been also successfully applied 
to the analysis of the edge profiles in the smectic phase formed by main-side chain 
liquid crystal polymers deposited on a solid substrate 58 ' 59 . In this experiment, the 
x-ray scattering has revealed the unsually large roughness at the surface 58 . Then 
the scanning tunneling microscopy has shown the clear sign of steps at the free sur- 
face, an indication of dislocations inside the sample 59 . The value of the parameters 
obtained in this experiments are as follows: 1 /y/KB = 5, A = sjKjB = 10lA« 5d. 
The smectic length A is much larger than the smectic period d, reflecting the fact 
that in the main-side chain liquid crystalline polymer we have two characteristic 
lengths. One is the size of the mesogenic unit related to d and one is the full length 
of the polymer chain related to A. We also see that K increases with the length of 
the polymer. 

So far we have assumed that steps at the free surface are due to dislocations 
inside the sample. However, it should be noted that steps at the surface do not have 
to be related to dislocations(Fig.22) and can have the same origin as steps in solids. 
Theoretically, the energy associated with a step 60 (Fig. 14) is an order of magnitude 
larger than the one associated with the dislocation (Fig. 7). The extra cost in the 
distortion energy per unit length shown in Fig. 22. is A'jd/2, where A7 is the surface 
tension difference for the PS and PBMA block of copolymer. It is roughly equal to 5d 
dyn/cm. At the same time the extra energy associated with dislocation is \J KBd/2 
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which is equal to 0.5d dyn/cm. An order of magnitude difference clearly shows that 
the configuration with dislocations (Fig. 7) is favored over the configuration with 
steps(Fig.22) in this particular system. Also TEM micrographs 22 clearly show that 
under the step at the surface we have indeed edge dislocations. However as we have 
discussed in Section II in some cases the configuration with steps can have lower 
energy than the configuration with dislocations. 

For small surface tension the edge profiles 49 strongly depends on K s . Fig. 23 
shows the results for 7 = 0. The most interesting feature of u eq (x, z = 0) is the 
nonmonotonic behavior of the profile as a function of x, observed for finite D. This 
is due to the finite size of the film and zero surface tension, since only in this case 
it is favorable to distort more the interface than the bulk. Large surface tension 
would make the profile perfectly monotonic. In very thin films we expect that the 
surface tension is important in determining the step structure when 7 a/ K s B / D. 
In thick films the surface tension dominates when 7 > \J~KB. We find that the 
width of the profile grows with increasing K s . 

VIII. Elastic interactions between dislocations 

So far we have considered the properties of a single dislocation inside the sys- 
tem. Now let us consider two edge dislocations in the infinite system: one located 
at (x, z) = (0, 0) and one at (x, z) = (/, h). Since the equations for the displacement 
u are linear the total displacement, u b due to dislocations 1 and 2 is the linear 
combination of individual displacements u h and u^K The distortion energy can 
be conveniently written (using Eq.(5.3)) in the following form: 

F b (h, I) =E™ + E™ + \B J dx dU ^ Z) (u£\x, z = 0~)- „W (x, z = 0+)) 

(8.1) 
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Inserting the solutions obtained in section V we find 10 ' 19 

F b {h, I) = + 4 2) ± ^EEl^I^!^ f dqexp (-\ q 2 h) cos (ql) (8.2) 

All J 

Here plus sign corresponds to like dislocations (same sign of the Burgers vector i.e. 
same sign in condition (5.6)) and minus sign to the unlike dislocations; i = 1, 2 
is the strength of the dislocations (see Eq(5.6)). The interaction energy is negligible 
for two edge dislocations located in the same plane (z = h = 0). Also for I 2 3> AXh 
the interactions are weak. Two like dislocations (same Burgers vector) always repel 
each other along the x direction (dF b /dl < 0) and attract each other along the z 
direction only for I 2 > 2Xh (dFb/dh > 0). For I 2 < 2Xh two like dislocations repel 
each other also along the z direction. 

For two screw dislocations located at x = and x = I one finds 43 that the 
interaction energy is zero in the harmonic approximation. If one includes anhar- 
monic terms 46 in a simple perturbation scheme one finds that two screw dislocations 
interact with a logarithmic potential (per unit length of dislocation line): 

F M = ±2 [B n -^ + K 3 (^ + ^)) .n( i/r „), (8.3) 

where r is the larger core size of two screw dislocations. According to Eq(8.3) 
two screw dislocations of same sign (plus sign in Eq(8.3)) attract each other with 
the long range potential. Two screw dislocations of opposite sign (minus sign in 
Eq.(8.3)) repel each other preventing their annihilation. Although the procedure 
applied to this problem 46 casts some doubts on the full validity of this equation, 
nonetheless it shows that the complete role of the anharmonic terms in the theory 
of dislocations is not known. 

So far we have considered the infinite system with no boundary effects. Here 
we go one step further and consider the finite symmetric system, bounded by two 
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free surfaces (Fig. 8a). The distortion energy is given by the bulk term (Eq.(5.4)) 
and two surface terms F s (see Eq.(6.4)). Minimizing this energy we find Eq.(5.3) 
and two boundary conditions for each free surface (see Eq(6.5)). The solution is the 
sum of individual solutions (since the equations are linear). Because in our system 
7 > \J~KB two edge dislocations are located in the same plane in the middle of the 
film and thus their bulk interactions are negligible. Now the only elastic interactions 
are surface induced. Inserting the equilibrium solution for two dislocations at a 
distance / apart we find the surface induced elastic interaction per unit length of 
dislocation: 

FS) = ~ ^2 J WW 1 ± e " ) ( 8 - 4 ) 
where the plus sign corresponds to like dislocations and minus to unlike dislocations. 

a x {q) 



/(«) = 



a 2 (<?)exp (q 2 XD) + a^q) 



and a\{q) and a,2(q) are given by Eqs.(6.11-12). The total distortion energy for 
these dislocations is given by the sum of F 8 (l) and Ff,(h = 0, Z) (Eq(8.4),Eq.(8.2)). 
For large a s (how large a s should be depends on the film thickness and the surface 
elastic constant, typically a s > 2 is sufficient) the interactions are repulsive for all 
distances for like dislocations and attractive for opposite dislocations. For a s 1 
we find from Eq.(8.4) the following form of the elastic interaction energy: 

Fs{l) ~ 16^ ( 6XP (8 ' 5) 



where £ = y/ a s DX/2 = a/ 7-D/2S is proportional to the size of the edge profile 
as could be expected. We have neglected the bulk contribution to the interaction 
energy since at large / it decays as exp ( — (l/d) 2 ). In the extreme limit of a s — > oo 
(which correspond to the system sandwiched between two walls) we find that this 
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energy diverges for like dislocations. For large a s we find for like dislocations: 

= a a VKBnUnVd? _ BnWn m dH , { ^ D) . (8 . 6) 
871-4 

and for unlike ones: 

F 8 (l) = Bn {1) n {2) d 2 l/(8nD). (8.7) 

As could be expected unlike dislocations attract each other while like dislocations 
repel each other in the system with large surface tension. 

For small ot s , but still larger than 1, the behavior of F s (l) is more complex, 
than in the previous case. For like dislocations the potential is repulsive at shorter 
distances and attractive at larger distances, which results in the minimum of F s (l) 
roughly at / ~ V\D ps dy/N (D = Nd). For a s = 1 we find 

F -<" * ( 2 ± ( 2 - -P<-' 2 /4AI») (8.8) 

This behavior of the interaction potential can be understood by comparing the 
edge profiles for large and small surface tension. In the former case the profile is a 
monotonic function of the distance (Fig. 20) and like dislocations repel each other 
minimizing the free energy associated with the surface tension. In the case of small 
surface tension, surface energy associated with curvature dominates, the profiles 
become nonmonotonic functions (Fig. 23) of the distance and like dislocations attract 
each other at large separations. For a s C 1 the attraction is even more pronounced. 

Finally let us discuss the interactions between two screw dislocations in a finite 
system in a harmonic approximation. The solution given by Eq.(5.14) satisfies 
also the boundary condition (6.5) at the free surface. Thus we can use it for the 
calculation of the distortion energy in the system bounded by two free surfaces 

(Eq.(6.4)) The surface induced distortion energy per unit length is as follows: 

F a (l) = 1 {n^) 2 d 2 /{nD) ln(L p /r^)+ 7 (n^) 2 d 2 /(wD) In (L p /rf ) 

(8.9) 

±2 1 n^n {2) d 2 /(nD)ln(L p /l) 
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Here L p is the horizontal size of the system and D is the thickness of the film. Bulk 
terms (5. 15), (5. 16) and (8.3) are not included in Eq(8.9). The like screw dislocations 
repel each other and the unlike dislocations attract each other, contrary to the case 
of bulk induced interactions (8.3). In both cases the potential is long ranged and 
has the same logarithmic form. The bulk term (8.3) is much larger than the surface 
term even for very small D ~ lOd. However it should be noted that (8.3) has 
been obtained from the anharmonic terms. At the harmonic approximation the 
bulk interactions are zero and then (8.9) is the only contribution to the interaction 
potential. 

IX. Fluctuations induced interactions between dislocations 

Dislocations are low-dimensional manifolds just like strings 61 , membranes in 
solutions 8 ' 62 , steps at the solid- vacuum interfaces 63 , polymers 64 , interfaces 65 and 
domain walls in monolayers 66 . This fact has not been fully recognized so far( see 
however 7 ' 67 ' 68 ). Here we would like to explore the consequences of the flexibility of 
dislocations on their behavior in finite lamellar films. Such studies are important 
for understanding the behavior of walls formed by dislocations in liquid crystals, 
diblock copolymers or lamellar phases of microemulsions, formation of meniscus in 
lamellar phases, interactions and mobility in systems of dislocations. 

From the seminal work of Helfrich 62 it follows that flexible strongly fluctuating 
objects such as membranes experience long range repulsive forces when brought 
together 8 . They are called undulation forces and arise from strong fluctuations and 
short range repulsion between the objects. In the case of membranes in aqueous 
solution this short range repulsion is due to the hydration forces 69 . Here we shall 
provide the same analysis and discuss these forces for edge dislocations. 

The single edge dislocation of length L y located along the y axis and fluctuating 
in the x — y plane (Fig. 24) in the bulk lamellar system can be described by the 
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following elastic energy: 

where u x is the displacement along the x direction (Fig. 24), and Eq is the line 
tension given by Eq(5.7). Due to strong fluctuations < u 2 x > grows as the length 
of the dislocation L y . However these fluctuations are damped in the presence of 
other dislocations. In a system of parallel dislocations located in the x — y plane 
and separated from one another by a distance /, the full elastic energy: 



H = \fdxdy\^f- 



where the dislocation compressional modulus B = finksT) 2 / (16E 1 3 ) arises self 
consistently from the fluctuation free energy, in the same way as in a stack of 
membranes 62 . Here / is the average separation of edge dislocations. Following 
Helfrich 62 we find, using Eq(9.2), the undulation interactions per unit length, be- 
tween the dislocations: 

Eq(9.3) differs from Helfrich result only by a numerical factor. As we discussed 
above, short range repulsion must be present in the system in order for the Helfrich 
mechanism to hold. In general the source of the short range repulsion can come from 
the core structure. The edge dislocation with the core splitted into two disclinations 
(Fig. 4c) is indeed characterized by large core energy and complex core structure and 
consequently two such dislocations of opposite or same Burgers vectors might repel 
each other at short distances. 

The fluctuations of the screw dislocations are more subtle than those of edge 
dislocations. In particular the line tension is not equal to the energy of a straight 

37 



screw dislocation per unit length (compare Eq(5. 19-20)), but is much larger 44 . The 
crucial point here is that in the process of deformation the screw dislocation aquires 
an edge character. This is the main source of the large configurational tension given 
by Eq(5.20). 

X. Unbinding transition for edge dislocations 

Our strongly fluctuating objects can also attract each other (e.g. membranes 
interact via the van der Waals potential) and at large distances this interaction 
may compete with undulation forces. At sufficiently low temperatures or suffi- 
ciently strong attraction the fluctuating objects stay close minimizing their energy. 
When the temperature is raised or attraction weakened they separate, maximizing 
their entropy. The transition from the bound to the free state is called the unbind- 
ing transition 70 . As follows, there are three elements necessary to induce unbinding 
transition: strong fluctuations of objects, short range repulsive forces and suffi- 
ciently strong attraction between the objects at large distances. We have already 
discussed them in the previous sections. Here we will show that edge dislocations 
can undergo the unbinding transition at suitable conditions 71 . 

Since for opposite dislocations the elastic interaction (Eq(8.7)) is attractive, 
this competes with the Helfrich repulsion given by Eq(9.3). In the extreme limit 
of a s — > oo (which correspond to the system sandwiched between two walls) we 
find that, as a result of this competition, the opposite dislocations are stabilized 
at the distance, l eq = Tr((3(kBT) 2 D)/(2d 2 BE )) 1 / 3 . This distance corresponds 
to the minimum of the sum of undulation and elastic interactions F H (l) + F s (l) 
(here for simplicity we assumed n^ 1 -* = = 1). One finds for the typical values 
of the smectic parameters and D = lOOjum that this distance is very small, i.e. 
l eq = lOOA. It is understandable that we cannot expect an unbinding transition in 
this case, since in the limit of hard walls the elastic interations grow linearly as a 
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function of /. For finite a s , although in principle the unbinding transition would be 
possible in thin films, we find for typical smectic liquid crystal parameters (7 = 30 
dyn/cm, \JKB = 6 dyn/cm, UbT = 4 • 10 -14 erg (room temperature), A = d = 
30A) that it does not occur. For all film thicknesses the opposite dislocations are 
stabilized at finite separations, providing that the Helfrich mechanism holds in their 
case. In general these results tell us that in the case of large surface tension elastic 
interactions are much stronger than the Helfrich repulsion. 

For small surface tension (comparable to V KB) the undulation interactions 
(Eq.(9.3)) compete with the attractive potential for like dislocations (Eq(8.8)) in 
thick films (large D). As the thickness of the film is increased the undulation 
interactions win and the unbinding transition occurs. 

In our calculations we have used the same parameters as above, assuming 
a s « 1, 2 > K s /(Kd) > 1 and rS 1 ^ = = 1. We have also assumed r c = d and 
E c = VKBd 2 /(2r c ). 

We have noted that practically for a s > a* = 1.5 no bound states for disloca- 
tions exist for any film thicknesses (for K s /K = 1). This critical value of a* depends 
on the surface bending elastic constant e.g. for K s /Kd = 10 we find a* = 4. 

In Fig. 25 the plot of F(l) = F s (l) + F H (l) versus / is shown for three film 
thicknesses (N =30,40 and 54) (D = Nd). As the thickness of the film is increased 
the value of the energy at the minimum at Iq ~ 3dVW approaches zero. The 
unbinding transition for K s /Kd = 1 and a s = 1 takes place at = 54, when 
F(Iq) = F(oo). For thicker films (N > 54) the dislocations are free. For thinner 
films (N < 54) they are bound. 

In Fig. 26 we show the diagram of K s /Kd versus N. The dotted line separates 
the regions where the unbound state is expected from the one where the dislocations 
are bound. For large surface bending elastic constant we find the unbinding tran- 
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sitions at larger thicknesses of the film. On the contrary for larger surface tensions 
we find the unbinding transition takes place for thinner films (Fig. 27). 

One notes that the transition is extremely sensitive to the value of the surface 
tension for 7 \JKB. Changing it by 1% we get the decrease of N by a factor 
of 2. We can check that at the room temperature changing the temperature by 
one degree upwards decreases the thickness at which the unbinding transition takes 
place (for a s = 1,K S /Kd = 1) from N = 54 to N = 52. Therefore there are four 
parameters which affect the unbinding transition: surface tension, surface bending 
elastic constant, film thickness and temperature. All of them can be controlled, 
and so our results can be tested experimentally, especially that the number of soft 
condensed matter systems characterized by small surface tension is quite large. 
XI. Unbinding transition for dislocation loops 

Dislocation loops in bulk smectic systems have been first studied by Kleman 72 . 
Later Helfrich 67 argued that the nematic-smectic (NA) phase transition in liquid 
crystals could be initiated by the unbinding of dislocation loops. The mechanism 
would be similar to the one proposed by Kosterlitz and Thouless (KT) for the phase 
transitions in two dimensions 73 . As shown by Nelson and Toner 68 the smectic with 
unbound loops behaves like a nematic. The estimate of the transition temperature 
in the loop mechanism is simple and similar to the KT estimate. The free energy 
of the loop in the bulk has the following form: 

F loop = F C L - k B T^ lnp, (11.1) 

where L is the total length of the loop, F c is the self energy of the loop per unit 
length and p is the number larger than one (for the problem considered on the 
lattice this number would be close to the coordination number of the lattice) . Here 
fce(lnp)/(i is the entropy of the loop per unit length; we count the number of loop 
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configurations and neglect the entropy of placing and closing the loop. The latter 
is proportional to the log of the size of the loop and in general controls only the 
density of loops but does not influence the transition temperature, 

T NA = F c d/{k B lnp). (11.2) 

At this temperature Fi oop changes sign and the spontaneous growth of loops occurs. 
Then, according to Nelson and Toner, the layered smectic structure is destroyed. 
A very instructive approach to this problem can be found in Ref.7. Here we shall 
estimate the influence of the boundaries on the loops, their growth and size. 

Let us consider a smectic liquid crystal film of thickness D, bounded by two 
surfaces located at z = ±D/2, with a single circular dislocation loop of radius R 
located in the middle of the film at z = 0. The loop is described by the following 
condition for the vertical displacement , u(r±,z), of smectic layers: 

u (t ± ,z = 0) = I ' , !!! r± !-o ; (11-3) 

v - L ' ' \ sgn(z)d/2, if|rx|<i2. v ' 

Since in the smectic liquid crystal the surface tension is larger than \JK~B (a s > 1) 
the dislocation loop is stabilized in the middle of the film. The solution of Eq(5.3) 
consistent with the condition (11.3) and satisfying the boundary condition (6.5) 
(and similar condition for the second surface obtained from (6.5) by changing the 
sign of z). is given by w e(J (r_ L , z), 

u eq (r ± , z) = u b (r ± , z) + u p (r_ L , z) (11.4) 

where 72 

u 6 (r ± ,z) = 8gn(z)^ / exp {-\q 2 \z\) 6XP J 1 (qR) (11.5) 

An J q 



41 



and 74 

u p (r_L,z) = ^ J cfcu 6XP ^ ±r± ^ J l (qR) (exp (Xq 2 z) - exp (-Xq 2 z)) f(q). 

(11.6) 

Here g = |qj_| and /(g) is denned after Eq(8.4). Inserting our solution into the 
distortion energy + F s (Eq(5.4,6.4)) we find the distortion energy of the smectic 
film of thickness D with the circular dislocation loop of radius R in the middle of 
the film: 

F(D, R) =2ttRE c + VKBd 2 -R 2 / dqqj\{qR) 

2 Jo . (11.7) 

POO \ ' 

-2VKB7id 2 R 2 / dqqjj(qR)f(q) 
Jo 

In the limit of 3> d we find the following form of the distortion energy F C (D) per 
unit length as a function of the film thickness, D: 

F C (D) = lim = Ec + _ VKSrf 2 r dg/( ? ). (11.8) 

The first two terms of the energy are associated with the bulk distortions (Eq(5.7)) 
and do not depend on the surface tension or the film thickness. In fact they are 
the same as for the edge dislocation. Thus the local distortions induced by a large 
dislocation loop are the same as the ones induced by an edge dislocation. The sum 
of E c and VKBd 2 /2r c is equal to F c (Eqs(ll.l,2)). The last term in Eq(11.8) is 
associated with the finite size of the system and the surface tension. In the limit of 
D — > oo it vanishes and consequently F C (D) approaches F c as it should. In order to 
estimate the NA transition temperature as a function of the film thickness we use 
F C (D) in Eqs(ll.l,2). In general, p should also change with D. For thin films this 
number should be close to 4 (coordination number for the two dimensional square 
lattice) while for thick ones it will be close to 6. Here we neglect this contribution 
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to Eq(ll.l). From Eqs(ll. 1,11.8) we find that 48 ' 74 : 



T NA (D)-T NA (oo) \fKBd 2 Kl(Kn ^ fllQl 

Tf, — 7— x = —7=—C{a a ,K a /{KD)). (11.9) 

Tna{oo) VXDF c 

where 

C(a s , rj) = - t°° dq- (1 ~ ~ ^g 2 ) , , 

and Tjva(oo) is given by Eq(11.2). The function C(a s , K S /KD), depends on 
K S /{KD) very weakly, thus the dominant dependence of Tna(D) on D is 1/y/D. 
For ck s = 5 (typical for smectic liquid crystal 51 ) we find C=1.25. For the purpose 
of rough estimate we neglect the core energy E c and set C = 1 in Eq(11.9), finding: 

T NA (D)-T NA (oo) = J2r^_ 

T N a(oo) v^d' 1 ' ' 

The transition temperature is larger in thin films than in the bulk, thus the smectic 
phase is stabilized in films and near free surfaces. The large surface tension and finite 
size effects are responsible for this behavior. For three layer film, D = 3d, A as d = 
30A, and the core radius r c equal to the width of the molecule, r c «i3A(chosing the 
smallest length scale in the problem is consistent with E c = 0; usually it is assumed 
that r c ~ d; however we may encounter situations when the size of the core can 
be smaller 39 or larger (see Fig.4b)), the transition temperature can be shifted by 
tens of degrees upwards in comparison to the bulk. It means that it should be 
relatively easy to overheat a smectic film. It must be realized however, that for 
large temperature shifts it is not justified to neglect the temperature dependence of 
the smectic elastic constants. 

Growth of the dislocation loops is only one particular mechanism of the nematic 
smectic phase transition. However it can be preempted by some other mechanisms 7 . 
In thin films different mechanisms can be characterized by different temperature 
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dependence on D. For the loop unbinding mechanism we have already obtained 
5T NA = (T NA (D) - T NA (oo))/T NA (oo) ~ 1/y/D. For the 3D XY universality 
class 75 ST NA ~ D~ 1 / v with 76 v ~ 0.7, by comparing the correlation length to 
the thickness of the film. Finally for the first order phase transition we find 77 
STna ~ l/D. Since for each model of NA transition we find different behavior of 
T NA {D), this temperature can serve as yet another test for the order and mechanism 
of NA transition. 

So far we have tacitly assumed that 7 is finite. If we set 7 — > 00, the surfaces 
become rigid and we obtain the film sandwiched between solid boundaries. The 
distortion energy of the loop located between solid boundaries is proportional to 
the area of the loop, contrary to the previous case of the freely suspended film when 
the loop distortion energy has been proportional to the length of the loop. Since 
F{D,R) grows now as R 2 while the entropy still as R, we conclude (see Eq(ll.l)) 
that the loop mechanism is supressed in this case for any finite D. Now let us 
estimate the size of the loop in the film contained between solid boundaries. At 
T > Tn A (oo) the loop can grow up to the point when Fi oop changes sign from 
negative to positive. This happens for the finite radius of the loop, R eq . For 
D,R^> d and 7 — > 00 we can make the following approximations: 



f 

Jo 



dqqJi(qR)f(q) » (11.12) 



and 

' dqqjKqR)^^-. (11.13) 

Here A ps 0.5 is a constant. Now using Eq(ll.l) and Eqs(11.7,8) we find that 
Fioop = at R = and at 

R - XDFc T - T na(oo) (uu) 

eq AVKBd 2 T NA (oo) y ' } 
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As we see for D — > oo we get R eq — > oo. Certainly before _R eg could grow to infinity 
the smectic phase would undergo a transition to the nematic phase. The mechanism 
of this phase transition would be different from the loop unbinding mechanism. 
However in the overheated smectic between solid boundaries one should observe 
the loops of size given by Eq(11.14). 

Our approach to the loop deformations neglected completely the fact that the 
loop can have an edge as well as screw parts. Certainly for large thicknesses of 
the film the screw character of the loops should be visible 78 . In the experiments 
conducted on lamellar phases of non-ionic surfactants the density of loops increased 
at the approach to the isotropic-lamellar phase transition 78 . In fact large number 
of dislocation loops and their interactions might modify the transition temperature. 
This is not taken into account in our approach. 

XII. Dynamical properties of dislocations: mobility 

Dislocations can move inside a sample under the action of an external stress or 
stress created by another dislocation. Let a E be this stress (excluding the stress 
created by the dislocation itself). The force F E acting on the dislocation (per unit 
length) can be calculated similarly as in solids from the Peach Koehler formula 33 : 

Ff = tiki(yf m b m tk (12.1) 

where b is the Burgers vector of the dislocation (of length |b| = nd) and t is the unit 
vector along the dislocation line. The stress tensor is given by Eqs. (4.6-8) (see also 
Ref[43]). This force sets the dislocation in motion. If the dislocation moves slowly 
with respect to the sound velocity the effective inertial forces can be neglected, since 
the effective mass is very small in this case (typically a row of molecules). Thus the 
equation of motion can be simplified to the following form: 

F E + F v = 0, (12.2) 
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where is the viscous force (section IV). Now based on Eq(12.2) we can discuss 
the motion of screw and edge dislocations. Since the distiction between these two 
types of dislocations are more pronouced in smectics than in solids we shall consider 
their mobility separately. Let us first consider the edge dislocations. There are two 
basic kinds of motion: the glide and climb. The former is perpendicular, while the 
latter is parallel to the layers. These two kinds of motion are different: glide is 
conservative, while climb requires transport of matter by diffusion (permeation). 
In solids glide is favored and in fact is responsible for plastic deformations at low 
temperatures (T < T m /2; T m is the melting temperature). By contrast climb seems 
to be favored in smectics (at least in thermotropic smectics) and can be induced 
by the compression or dilation of layers while glide requires bending of the layers. 
Let us first consider the climb motion of an edge dislocation in a smectic system 
submitted to the normal stress af z = a (Fig. 28). The force acting on the dislocation 
of Burgers vector length |b| = nd is (Eq(12.1)): 

F E = -and, (12.3) 

while the friction force is 

F v 

From Eq(12.2) we find the velocity: 

v = —ma, (12-5) 

where m is the mobility of the dislocation, a priori being a function of the velocity, 
the length of the Burgers vector, nd, and the material constants. The mobility of 
an elementary dislocation 79 is related to the self diffusion coefficient, D\\, for ther- 
motropic smectics assuming that the mean distance between jogs, I, along the line 
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ndv 



m 



(12.4) 



is comparable to the molecular size. This assumption is reasonable in thermotropic 
smectics and gives with analogy to solids (metallurgical model): 

m ~^rT (12 - 6) 

Here v mo \ is the molecular volume and / ~ d. This model is not correct for giant 
dislocations, especially those with core splitted into two ±1/2 disclinations. In 
this case it is better to assume that the dislocation acts as an obstacle of width nd 
perpendicular to layers (Fig. 29). The flow of a smectic around such a ribbon leads to 
four permeation boundary layers in which dissipation occurs due to permeation 80 ' 81 . 
A calculation of dissipation gives the force of friction acting on the ribbon and the 
following expression for the mobility of the dislocation: 




m = W— , (12-7) 



where A p is the permeation coefficient and \i is the shear viscosity (see section IV) . 
A more sophisticated calculation 82 taking into account deformations of the layers 



and assuming strong permeation (l p = a/A p/ u > A = y/K/B) gives the following 
formula: 

m =4fe (i2 - 8) 

where A is a numeric factor of the order of one. The mobility does not depend on 
the Burgers vector nor on the velocity, due to the fluidity of the layers. The two 
formulas (12.6) and (12.8) are compatible since A p ~ D\\v mo \/kBT (section IV) and 
\i ~ k B T/W\\. 

Let us now describe the glide motion of an edge dislocation. In this motion, 
the dislocation moves perpendicularly to the layers (Fig. 30). Two cases must be 
considered. Below the Peierls-Nabarro stress the dislocations remain pinned in its 
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Peierls valley (at zero temperature) . For smectics this stress has been calculated by 
Lejcek 83 : 

apN = -B^~^— exp (—iir c /nd). (12.9) 

The process of glide must be thermally activated for af z < &pn- It requires the for- 
mation of small bulges (Fig. 31) which can afterwards spread sideways. For stresses 
larger than <jpn dislocations can glide easily and the mobility is 84 



16r cV /7rAr e 

771 ~ vnd(l + 3\/(2r c ))- [ UJ 

The glide mobility is of the same order of magnitude as the climb mobility Eq.(12.8). 
The more dissociated the core is (r c large, Fig. 4b) the easier glide proceeds (large 
mobility). However for the giant dislocations with core splitted into ±1/2 disclina- 
tions the glide is much more difficult than climb, since the core is extended along 
the x direction. 

To finish this section let us give the mobility of screw dislocations. These 
defects can easily glide without the layer resistance. The mobility is given by the 
following formula: 

8n 2 r 2 

m = —-±. (12.11) 

find 

The motion of the screw dislocation is conservative and thus permeation is negligi- 
ble. The friction force is due to the viscosity of the medium \i. 

XIII. Microplasticity and helical instability of screw dislocations 

Experimentally we can reveal the dynamical properties of dislocations by com- 
pressing a homeotropic wedge sample (see section II) and recording its viscoelastic 
response. More precisely one can impose a step like variation or a sinusoidal varia- 
tion of the thickness with piezoelectric ceramics and measure the normal stress as a 
function of time. The first experiment of this type was performed by Bartolino and 

48 



Durand 85 in smectic-A liquid crystal. By applying very small thickness variation 
(a few A or less) they showed for the first time the elastic behavior of smectics 
in compression normal to layers and its plastic relaxation in time. Ribotta also 
studied with this system the undulation instabilities of layers in dilation and their 
subsequent instabilities leading to focal parabolas 34 ' 86 . In order to study disloca- 
tionsand plasticity a new cell with stackings of piezoelectric ceramics, allowing for 
much larger deformations, was used 87 . The step like thickness variations (ranging 
from d to lOd typically) and the normal stress as a function of time are shown 
schematically in Fig.32ab. The stress first jumps abruptly to ao ~ B5o/d (elastic 
regime) and then relaxes exponentially to zero with the characteristic time r. This 
time depends on the initial stress 88 (Fig. 32c) and decreases by successive jumps 
with the increase of cxo- Each jump is associated with the change of the thickness 
by N p = pNi layers (p is the index of the jump) and by the relaxation time r p 
defined by: 

111 , 

- = - + -. (13.1) 

The time r g is the relaxation time at small deformations. The experiment shows 
that Ni ~ 3 independent of the thickness of the sample D and the angle a between 
the plates. One finds from the experimental results that r g ~ D/a and t p ~ D/p, 
where p is the index of the jump (Fig. 32c). Also t p is independent of a. Two 
relaxation mechanisms are involved in the observed phenomena. The first one is 
the motion of edge dislocations which are present in the sample because of the 
wedge geometry. The relaxation time for this motion is 

D 

Tq = —— 13.2 

9 mBa v ' 

proportional to D/a in agreement with experiment. From (13.2) we can extract 
the mobility. For 8CB at room temperature one finds m ~ 5 • 10 _7 cm 2 s/g and 
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from Eq.(12.7) we calculate A p 5 • 10~ 13 cm 2 /poise {ji 2poise). This gives 
/ p = a/7^V ~ 10A. This length is of the order of molecular size as expected from 
Eq.(12.6) and Einstein relation for viscosity. 

The second mechanism which explains the jumps of the relaxation time is the 
sudden instability of screw dislocations under the applied deformation 44 . If 

AH ± screw , . 

a >p — (13.3) 

then, each elementary screw dislocation (n = 1) joining two plates (and anchored 
to them) is unstable with respect to the formation of a helix of pitch D/p. This 
helical shape is equivalent to the removal (under compression) of p layers inside 
the cylinder in which the screw dislocation is located. (Fig. 33). That is why the 
relaxation time is inversely proportional to the jump order p. Finally we can use 
the formula for iVi: 

Nl = 2tt^^ (13.4) 

to estimate E c using the theoretical value of T screw (5.19-20). One finds E c ~ 
O.lBd 2 . This value of the core energy is compatible with the one obtained for the 
core filled with nematic. The relaxation time associated to jump p can also be 
calculated: 

(13.5) 



p 2n 2 mT screw p 

where m is the screw dislocation mobility given by Eq.(12.11). From Eq(13.5) and 
t p measurements we find m ~ 10 _6 cm 2 s/g in agreement with Eq(12.11). 

So far we have assumed that both relaxation processes are independent, but 
we know that edge dislocation must cross the screw dislocation in the process. This 
crossing leads to the formation of the jogs on the edge dislocation (Fig. 34). The jog 
has the screw character so its energy (Eq(5.15)) is 

*W ~ |J (13.6) 
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assuming n = 1 for both dislocations and r c = X = d. Experimentally B ~ 
5 • 10 7 dyn/cm 2 and d ~ 30A, thus E ]Qg ~ 10~ 2 k bT. This jog can be very easily 
thermally activated thus elementary dislocations can very easily cross. Also two 
edge dislocations interact only when they are in the same slip plane which is rarely 
the case, consequently each dislocation moves independently from the other dislo- 
cations, contrary to the case of dislocations in solids. Of course once again fluidity 
of the layers is responsible for this behavior of dislocations. Finally, one may ask 
why helical instability of screw dislocation is more favorable than the nucleation 
ex-nihilo of dislocation loops. The answer is as follows. The length increase of 
the screw dislocation when it aquires a helical shape of radius R is 2ttR x nR/D, 
whereas the nucleation of the dislocation loop of equivalent radius requires to make 
a length 2ttR of edge dislocation. The helical instability is thus much more favorable 
energetically than the nucleation ex-nihilo of dislocation loop. 

Dislocations may also play an important role in the process of growth. This 
has been observed at the smectic A-smectic B interface. When the smectic B grows 
from the undercooled smectic A, it generates stresses that are mainly relaxed by 
the motion of dislocations in the smectic A phase. This phenomenon is detectable 
owing to the entrapment of dust particles which are dragged away when they climb 
parallel to the layers. If the growth rate is too large, then dislocations may no 
longer relax the stresses efficiently and then smectic A breaks with the formation 
of many focal conic domains 84 . 

So far we have considered measurements performed on smectic liquid crystals. 
Similar measurements have been done for self assembling amphiphilic systems of 
C12E5 (pentaethylene glycoldodecyl ether) and water. Surprisingly this lamellar 
system orients very easily and spontaneously in homeotropic anchoring. For exam- 
ple few minutes are enough to obtain a good homeotropic sample without visible 
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defects, whereas several weeks of annealing are necessary to obtain similar results 
with phospholipid lamellar phase 26 . This simple observation shows that permeation 
is much faster in C12E5 than in usual lyotropic systems, the fact also confirmed by 
the direct measurements of the edge dislocation mobility 89 . The mobility has been 
found of the same order of magnitude as in thermotropic liquid crystals considered 
previously in this section. The screw dislocations are very numerous in this ly- 
otropic system 90 and presumably they favor permeation. In fact it was shown that 
these defects act as vortices when a pressure gradient is imposed perpendicularly 
to the layers. This mechanism leads to the effective permeation coefficient 91 

where L is the mean distance between dislocations and 

A = 0.07 + 0.009 In L/r c . (13.8) 

Finally we note that dislocations strongly affect flows as will be shown in the 
next section. 

XIV. Lubrication 92 

It is well known that lamellar phases are very good lubricants. As an example 
one can cite soaps, graphite, and also lamination oils, which under high pressure 
and temperature acquire a lamellar structure. In the classical theory, lubrication 
occurs when a fluid is sheared between two planar surfaces making a small angle 
a. The lubrication force is due to the increase of pressure in the locally compressed 
fluid. By replacing the isotropic fluid with smectic A one expects similar, but more 
spectacular effects, since apart from the easy flow of the layers past each other, there 
is a genuine solid like strength perpendicular to layers. Fig. 35 shows the lubrication 
geometry. The sample is sandwiched between two plates with homeotropic anchor- 
ing at an angle a. In this sample as we know there is an array of dislocations at a 
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mean distance d/a. At time t = the upper slide starts to move with the horizontal 
velocity v. At time t the thickness of the sample has varied perpendicularly to a 
fixed position by an amount AD = — vta, which is positive (dilation) when v < 
and negative (compression) when v > 0. This thickness variation creates a normal 
stress af z which makes the edge dislocations to climb together in a direction which 
favors the relaxation of the imposed deformation. The dislocations are dragged by 
the flow with velocity v/2 at time t = and then progressively slow down until they 
stop in the stationary regime (t — > oo). A direct calculation of the normal stress 
gives: 

a ^ = "2^ (1 " eXp( " t/T " 3)) ' 
with T g given by (13.2). The stress exerts a normal force on the slides, which adds 
to the one exerted by the hydrostatic pressure. Both terms can be calculated in 
the limit of a very small angle a. In the limit t ^> r g one finds 92 the force per unit 
length 

1 / L \ 3 Lv 



F n = -^hydrostatic + -^smectic = ~H \Jj j VOL + (1^-2) 

Here L is the horizontal size of the sample. The first term is related to the hy- 
drostatic pressure and is always present, even for isotropic lubricants, whereas 
the second one is specific for layered structures. In fact the second term can 
be much larger than the first one. We estimate -Fhydrostatic = 10 3 v[dyn s/cm 2 ] 
and -Fsmectic = 10 6 i>[dyn s/cm 2 ], assuming L = 2cm, D = 200/um, \i = lpoise, 
a = 2 ■ 10~ 4 and m = 10~ 6 cm 2 s/g. Since the normal force is three orders of magni- 
tude larger in the case of smectics it means that smectics better keep the two moving 
surfaces apart and consequently are better lubricants than an ordinary fluid of the 
same viscosity. Of course the full comparison of the lubrication properties involves 
the calculation of the force of friction. Assuming simple flow pattern v x = v z/D we 
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find the force of friction (per unit length) 

_ vL aLv . . 

^friction = -fJ>-pr ~ -. , (14.3) 

D Am 

and therefore the smectic behaves as a fluid with apparent viscosity 

"™> = " + 4^7 * " + "17' (14 ' 4) 

Thus the corrections to viscosity \i can be very significant even for very small angles 
a, because D is much larger then the permeation length l p ~ d (see section XII 
before Eq.(12.8)). For instance /U app = when a = 2 • 10 _4 rad, D = 200/um and 
/ p = 25A. This larger viscosity is of course not favorable for lubrication. Another 
consequence of this calculation is that lubrication effects must be taken into account 
when one measures the shear viscosity with rheometer. Nevertheless, it is possible, 
at certain conditions, to completely eliminate the lubrication effects and to measure 
the intrinsic shear viscosity of the sample. Even in this limit, the viscosity appears 
to be strongly dependent on the shear rate S = v/D; large at small S, fi app ap- 
proaches n at large S (Fig. 36). Careful observations of the motion of dust particles 
in the samples show that the viscosity increase at smaller shear rates is related 
to strong departure from the linear velocity profile (Fig. 37). In order to explain 
the rigidification of the velocity profile in the middle of the sample a mechanism 
of strong interaction between screw dislocation and the flow has been proposed 93 . 
Their stationary shape can be obtained by comparing the viscous torque that the 
flow exerts on the line with the elastic torque due to the line tension. The disloca- 
tions bend while acquiring a mixed character and tend to occupy a finite fraction 
of the whole thickness of the sample, when the shear rate becomes larger than the 
critical shear 

_ lesrf 2 ^ 16K 

bc ~-^~w (14 - 5) 
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The larger the shear rate the thinner is the zone where the flow is rigidified (small 
velocity) (Fig. 37). At very large S a dislocation is completely stretched and do not 
longer contribute to the total dissipation. In this limit the apparent viscosity // a pp 
is equal to the intrinsic viscosity \i. This model allowed us to explain most of the 
experimental results and to estimate the density of screw dislocations (10 6 cm~ 2 ) in 
usual sample of thermotropic smectic liquid crystals 94 . 

Summary 

In this review paper we have tried to summarize the progress in the theory of 
dislocations which took place in the recent decade. Both equilibrium and nonequi- 
librium, bulk and surface properties have been discussed. It has been shown that 
theoretical and experimental studies are very closely related in this field and that 
many interesting aspects of the behavior of dislocations can be studied experimen- 
tally. Many properties of dislocations in smectics are very much different from 
those of solids. We have tried to emphasize all the differences. In general we can 
say that smectics have certain solid like and liquid like properties combined in the 
very unsual way and in this review we have shown how the solid-like and liquid-like 
behavior of smectics influences dislocations. 
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Figure captions 

Fig.l Schematic picture of smectic (lamellar) liquid crystal obtained in computer 

simulation of a hard rod system by D.Frenkel. 
Fig. 2 Lamellar phase in binary mixture of water and surfactant. The surfactant 

molecuels form bilayers. From F.B. Rosevear J. Soc. Cosmetic Chemist 19, 581 
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(1968). 

Fig. 3 Schematic picture of lamellar phase in diblock copolymer system. The linear 
copolymer consists of two homopolymers A (solid line) and B (dashed line) 
joined by chemical bond (dot). Alternating A rich and B rich domains are 
shown. After Ref.6. 

Fig. 4 (a) The edge dislocation with small core, (b) The edge dislocation with ex- 
tended core (c) The edge dislocation with core splitted into two disclinations 
(d) The screw dislocation, after Ref.16 

Fig. 5 (a) Schematic cross section of a sample containing a dislocation array, (b) 
strain versus position (c) accompanying variations of tilt angle versus position 
(d) Lagerwall subgrain boundary at the smectic A smectic C phase transition, 
after Ref.13. 

Fig. 6 Scan of the well annealed sample in the lamellar phase showing regular steps 
due to the edge dislocations. Each step corresponds to the termination of the 
single bilayer. The shot noise level is approximately 10% of the step size, after 
Ref.26. 

Fig. 7 (a) AFM image of the holes and islands at the free surface of AB diblock copoly- 
mer (here PS/PBMA diblock copolymer). The height of the circular elevation 
is 310 A. (b) Cross section of the island, (c) The schematic representation of the 
structure with the dislocation. Since the islands are circular here we have dis- 
location loops. Locally the distortions induced by large edge dislocation loop 
are the same as those induced by linear edge dislocation(Fig.4a). The dark 
and white regions corresponds (see also Fig. 3) to PS and PBMA domains. The 
dark lines represent PS/PBMA interfaces, while dotted lines are Active PS-PS 
or PBMA-PBMA separations, (d) TEM picture showing the edge dislocation, 
after Ref.21. 
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Fig. 8(a) The idealized smectic meniscus of the freely suspended smectic liquid crystal 

film, after Ref.23. 
Fig.8(b) The arch-texture, after Ref.23. 
Fig.8(c) Knots after Ref.23. 

Fig. 9 Profiles of 80CB droplets obtained by Michelson interferometry at different 

temperatures. The indicated temperatures are the numbers in of °C below 

T N a- The apparent droplet diameter is 177 ± 5/um . The matching ofthe facet 

with the curved part of the droplet is tangential. 
Fig. 10 Possible droplet configurations. In (a) the layers are parallel to the substrate. 

The top surface shows a single facet parallel to the substrate and steps. In 

(b) the top layer is curved to follow the free surface and dislocation loops run 

through the bulk of the smectic liquid crystal. 
Fig. 11 (a) Part of the "goutte a gradins" seen through themicroscope in reflection with 

natural light, (b) Schematic view of a "goutte a gradin" (c) Probable structure 

of a "gradin" . 

Fig. 12 Droplet profiles of 40.8 for different temperatures below T N a (droplet diameter 
145/um). When the temperature decreases the droplet shows a clear tendency 
to develop a discontinuity in slope at the facet edge. A secondary facet is well 
visible on the profile at T NA - T = 9.8°C. 

Fig. 13 TEM picture of screw dislocations in lecithin (courtesy of M.AUain). Arrows 
show the emergence points of screw dislocations. Sample has been frozen and 
then fractured under vacuum. A replica of its surface is observed by TEM. 

Fig. 14 The smectic order parameter e, nematic order parameter s and the local angle 
normalized by their bulk values, as a function of the distance from the 
center of the screw dislocation core. The distance is measured in the units of 
the correlation length for smectic This length is of the size of a molecule for 
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small temperatures but diverges at the approach to the smectic-nematic phase 
transition temperature (in the case of continuous transition in liquid crystals 7 ). 
After Ref.42. 

Fig. 15 The screw dislocation tilted with respect to the layers. 

Fig. 16 Two configurations of a dislocation pinned at the solid substrate (a) unfavorable 
(b) favorable. 

Fig. 17 The energy of a single edge dislocation F* = (F(h)-E )8nX s /(VKBd 2 ) (A s = 

versus the distance from the surface h* = 2hK/K s for 7 j\/KB = 
(dashed line) and 'j/y/KB = 0.5 (solid line). In the latter case the shallow 
minimum (not visible on the scale of the Figure) is located at h* = 2.55. 

Fig. 18 The equilibrium location of the dislocation inside the film h* eq = h eq /d versus 
the surface bending elastic constant K* = K s /Kd for D/d = 3 - long dashed 
line; D/d = 10 - short dashed line; and D/d = 00 - solid line. a s = 0. 

Fig. 19. Clustering of dislocations in a wedge shaped sample (Nallet and Prost 55 ). 

Fig. 20 The edge profile u* = —u eq (x,z = 0)/d, versus x* = (x — l)/\/~Xh at the free 
surface of thin lamellar film deposited on solid substrate (dashed line). Here 
K s = Kd : 7/ y/KB = 5, D = 3d, h = d/2. For comparison (solid line) we show 
the profile, obtained in the previous approach 19 with finite size and surface 
effects neglected (D — > 00, 7 = 0, K s = 0). 

Fig. 21 A plot of the local film thickness against position measured along an axis paral- 
lel to the substrate and normal to the domain edge (see Fig. 5) (a) the average 
film thickness, D is d < D < 2d; (b) 2d < D < 3d; (c) 3d < D < Ad; 
(d) 5d < D < 6d. Dots display experimental results and the solid line was 
obtained from the theory (however here K s = 0). After Ref.22. 

Fig. 22 The same schematic picture as shown in Fig. 7 with steps at the free surface, 
without dislocations in the bulk, (this configuration is in fact unfavorable 
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in comparison to film structure with steps induced by dislocations Fig. 5 see 
Section VI). After Ref.60. 

Fig. 23 The edge profile u* = —u(x, z = 0)/d versus x* = (x — l)/V\d for the film on 
a solid substrate bounded by the interface with the surface tension 7 = for 
three thicknesses of the lamellar film: D/d = 3 - long dashed line; D/d = 10 - 
short dashed line; and D/d = 00 - solid line, (a) K s = Kd (b) K s /Kd = 10. 

Fig. 24 The top view of the z = surface containing the edge dislocation. The dashed 
line represents the unperturbed position of the dislocation and u x (y) is the 
amplitude of thermal distortions of the dislocation in the x-y plane (see also 
Eq(9.1)). 

Fig. 25 The dimensionless interaction energy per unit length V(x) = 2n 2 \^N / (^KBd) 
(F s (x) + F H (x) -F a (0)) versus the separation x = l/VXdN for N = 30, A = 40 
and N = 54. Here N is the number of lamellar layers. As N increases the 
minimum of V(x) (at xq ~ 3) moves upward. At N = 54 V{xq) = V(x = 
00) = and the unbinding transition occurs (K s /Kd = 1, a s = 1). 

Fig. 26 The surface bending elastic constant K* = K s /Kd (a s = 1) versus the number 
of lamellar layers, N, at the unbinding transition. The dotted line divides the 
diagram into two regions. Above the line the dislocations are free and below 
they are bound. 

Fig. 27 The surface tension a s = ^/yKB (K s /Kd = 1) versus the number of lamellar 
layers, N at the unbinding transition. The dotted line divides the diagram into 
two regions. Above the line the dislocations are free, below they are bound. 
At the line the unbinding transition occurs. 

Fig. 28 The climb motion of an edge dislocation submitted to a normal stress af z . 
v is the velocity, b is the Burgers vector and t is the unit vector along the 
dislocation. 
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Fig. 29 (a) (b) climb of the giant dislocation and equivalent rheological model. 
Fig. 30 The glide motion of a dislocation submitted to shear (see Fig. 29). 
Fig. 31 Thermally activated bulge. 

Fig. 32 The schematic view of the dynamical measurements, (a) step like deformations. 

Initial thickness variation is 8q. (b) Normal stress measured experimentally. 

The characteristic relaxation time is r. (c) Time r versus cr m a 400/um thick 

sample of 8CB at room temperature. 
Fig. 33 Helical instability of screw dislocation. 

Fig. 34 Crossing of an edge and a screw dislocation with the formation of the jog. 
Fig. 35 Lubrication geometry. 

Fig. 36 The apparent (effective) viscosity versus the shear rate (8CB, D = 100/um, 
22°C from Ref.91). 

Fig. 37 Velocity profile in the sample and a dislocation distorted by the flow. 
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